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Thomas Schelling proposed an influential simple spatial model to illustrate 
how, even with relatively mild assumptions on each individual's nearest neigh- 
bor preferences, an integrated city would likely unravel to a segregated city, 
even if all individuals prefer integration. Aggregation relates to individuals 
coming together to form groups and global aggregation corresponds to segre- 
gation. Many authors assumed that the segregation which Schelling observed 
in simulations on very small cities persists for larger, realistic size cities. We 
devise new measures to quantify the segregation and unlock its dependence on 
city size, disparate neighbor comfortability threshold, and population density. 
We identify distinct scales of global aggregation, and show that the striking 
global aggregation ScheUing observed is strictly a small city phenomenon. We 
also discover several scaling laws for the aggregation measures. 
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1 Introduction 



In the 1970s, the eminent economic modeler Thomas Schelling proposed a simple space-time 
population model to illustrate how, even with relatively mild assumptions concerning every in- 
dividual's nearest neighbor preferences, an integrated city would likely unravel to a segregated 
city, even if all individuals prefer integration (|6l|7l[8l|9l). His agent based lattice model has be- 
come quite influential amongst social scientists, demographers, and economists, and a number 
of authors are testing the Schelling model using actual population data (|T3l l2ini;IT4l). The only 
quantitative analysis of such models we could locate in the literature are ([3) and (TV). 

Aggregation relates to individuals coming together to form groups or clusters according to 
race, and global aggregation corresponds to segregation. Many authors assume that the striking 
global aggregation observed in simulations on very small ideal "cities" persists for large, realis- 
tic size cities. A recent paper (ITOl) exhibits final states for a small number of model simulations 
of a large city, and some final states that do not exhibit significant global aggregation. However, 
quantification of this important phenomenon is lacking in the literature, presumably, in part, 
due to the large computational costs required to run simulations using existing algorithms. 

In this paper, we devise new measures to quantify the aggregation and unlock its dependence 
on city size, disparate neighbor comfortability threshold, and population density. We developed 
a highly efficient and fast algorithm that allows us to compute meaningful statistics of these 
measures(see Figures 3|8|10 ). We identify distinct scales of global aggregation, and we show 
that the striking global aggregation Schelling observed is strictly a small city phenomenon. We 
also discover several remarkable scaling laws for the aggregation measures. 



1.1 Description of the Model 

We expand Schelling's original model[^to a three parameter family of models. The phase space 
for these models is the N x N square lattice with periodic boundary conditions. We consider 
two distinct populations, that, in Schelling's words ® , refer to "membership in one of two ho- 
mogeneous groups: men or women, blacks and whites, French-speaking and English speaking, 
officers and enlisted men, students and faculty, surfers and swimmers, the well dressed and the 
poorly dressed, and any other dichotomy that is exhaustive and recognizable." We denote by 
B (black squares) and R (red squares) these two populations. See Figure 1. Together these 
agents fill up some of the N"^ sites, with V remaining vacant sites (white squares). Each agent 
has eight nearest neighbors (Moore neighborhood). Fix a disparate neighbor comfortability 
threshold T E {0,1, ... ,8}, and declare that a B or R is happy if T or more of its nearest eight 
neighbors are B's or R's, respectively. Else it is unhappy. 

Demographically, the parameter N controls the size of the city, v = V/N"^ controls the 
population density or the occupancy ratio (4), and T is an "agent comfortability index" that 
quantifies an agent's tolerance to living amongst disparate nearest neighbors. 

'Different authors frequently consider slightly different versions of Schelling's original model, i.e., different 
ways of moving boundary agents. All versions seem to exhibit the same qualitative behaviors, and thus we refer to 
the Schelling model. 
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We begin the evolution by choosing an initial configuration (described in Sect. 3) and ran- 
domly selecting an unhappy B and a vacant site surrounded by at least T nearest B neighbors. 
Provided this is possible, interchange the unhappy B with the vacant site, so that this B becomes 
happy. Then randomly select an unhappy R and a vacant site having at least T nearest neighbors 
of type R. Provided this is possible, interchange the unhappy R with the vacant site, so that R 
becomes happy. Repeat this iterative procedure, alternating between selecting an unhappy B 
and an unhappy R, until a final state is reached, where no interchange is possible that increases 
happiness. For some final states, some (and in some cases, many) agents may be unhappy, but 
there are no allowable switches. 



2 Schelling's results 



Schelling considered the case = 8,T = 3, and v = 33%. The details of Schelling's simu- 



lations (including the algorithm for choosing initial conditions) are presented in Sec. 6.1 in the 
supporting materials. The final state of a typical run of Schelling's original model system is 
presented in Figs. [T] Schelling performed many simulations by hand using an actual checker- 
board, and observed that the final states presented a significant degree of global aggregation. 
He equated the global aggregation with segregation of a city. 





Figure 1: A: A simulation of Schelling's original model with 
N = 100. 



8; B: Our simulation with 



In this paper, we investigate whether the global aggregation that Schelling observed for 
very small lattices persists for larger lattices. In Fig.[lp , we present a characteristic final state 
for our simulations with = 100. Comparing Figs. [l]\ and [ip , one can see a striking 
qualitative difference between the two final states. While there is some local aggregation in 
the final state with = 100, there is no global aggregation. By viewing the plots of this and 
other final states, one immediately sees that the global aggregation that Schelling observed is 
a small lattice phenomenon. In Sections 3 through 5 we quantify the global aggregation using 
several new measures, and we carefully analyze the structure of the final states for different 
values of T and v. In particular, we quantify the claim that global aggregation is a small lattice 
phenomenon. 
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3 Simulations 



We study the dynamics for large lattices and present our results for = 100. Figures 2-10 
(except Fig. |4]) are all based on = 100. In the supporting materials, we discuss the cases 

= 50 and N = 200, and claim that greater than 100 does not lead to qualitatively or 
quantitatively different states and phenomena. We restrict our discussion to cities having an 
equal number of B's and R's. A separate manuscript (.12) will study the dynamics with different 
proportions of B's and R's. 

We consider T = 3, 4, 5 and v between 2% and 33%. The system does not evolve very much 
for other values of T: for T = 1,2 almost all of the agents are satisfied in most of the initial 
configurations, while for T > 6 there are almost no legal switches. Values of v larger then 33% 
correspond to unrealistic environments. For each pair of parameters T and v, we perform 100 
simulations (statistics are presented below and in the supporting materials). This number of 
simulations was chosen to ensure a 95% confidence interval for parameter estimation. 

We choose the initial configuration by starting with a checkerboard with periodic boundary 
conditions. Demographically, a checkerboard configuration is a maximally integrated config- 
uration. We then randomly remove V/2 of both B's and R's (thus keeping equal numbers of 
both agents). We permute agents in two 3x3 blocks. Alternatively, we could choose a random 
initial configuration. In general, except for small v, the final states are quantitatively similar as 
for our Schelling initial conditions. 

In Fig. [2} we present characteristic final states for different values of T and v. Visually, the 
aggregation in the final states along each column (with fixed v) is substantially different than in 
the final states along each row (with fixed T). More characteristic final states are presented in 
Figs. |5]j7]in the supporting materials . We begin the next section by defining several measures 
of aggregation that enable us to quantify this observation. 

4 Analysis 

4.1 Measures of Aggregation 

One can see from Fig. |2] that aggregation is a multifaceted phenomenon, and thus requires 
several measures to describe. Schelling used two measures to quantify the aggregation: (1) the 
average over all agents of the quotient of the number of like to unlike neighbors, and (2) the 
number of agents whose eight nearest neighbors all have the same label, and thus are completely 
surrounded by like agents. We call the latter quantity seclusiveness and denote by A'q. 

We introduce three additional measures. Aggregation manifests itself in two ways: (i) re- 
duced number of contacts between agents of different kinds and (ii) the dominance of agents of 
either kind in a certain area. The first of the following measures of aggregation quantify (i) and 
the second and third measures quantify (ii). 

(3) The adjusted perimeter per agent p of the interface between the different agents suitably 
adjusted for the vacant spaces. The adjusted perimeter is defined as twice the total number of R- 
B connections plus the total number of connections between R and B agents with vacant spaces. 



4 




Figure 2: Characteristic final states for different values of T and t>: A: T = 3, f = 2%, B: 

T = 3,v = 15%, C:T = 3,v = 33%, D: T = A, v = 2%, E: T = A, v = 15%, F: T = 4, 
f = 33%, G:T = 5,v = 2%, H: T = 5, = 15%, l:T = 5,v = 33%. 
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all divided by the total number of agents N"^ (see the supporting materials). Demographically, 
p is the average number of contacts an agent has with the opposite kind or with vacant sites. 

(4) The scale, or maximum diameter of the connected components (which we henceforth call 
clusters) L. The measure L is defined as the side length of the smallest square needed to cover 
every cluster. For states consisting of compact clusters, larger values of L correspond to larger 
scales of aggregation. Checkerboard configurations and configurations consisting of compact 
clusters are two extremes; for the former L = N. 

(5) The total number of clusters in a configuration Nc- This intuitively appealing measure of 
aggregation is useful to describe final states having mostly large compact clusters. To see its 
limitation, observe that "the maximally integrated" checkerboard configuration with v = Q has 
just two clusters. 

The definition of p was motivated by analogies of these models with the physics of foams. A 
key observation is thatp is a Lyapunov function, i.e., a function defined on every configuration, 
and which is strictly decreasing along the evolution of the system until it reaches a final state (see 
proof in the supporting materials). Thus the system evolves to minimize the adjusted perimeter 
between the interface of the R and B agents. The final states are precisely the local minimizers 
of the Lyapunov function, subject to the threshold constraint. This Lyapunov function is also 
the Hamiltonian for a related spin lattice system related to the Ising model ([5]). 

In Fig. |3| we plot average values of aggregation measures (2)-(5) introduced above for the 
final states with T = 3, 4, 5 and several values of v. The linear relationships of these disparate 
aggregation measures on population density seems remarkable. 

4.2 Global Aggregation Dependence on T 

From Fig. |2} one observes that: (i) the final states with T = 3 are very sparse, with a great deal 
of interweaving between both kinds of agents and vacant spots; (ii) the final states with T = 4 
consist of compact clusters (that look like "solid" objects); and, finally, (iii) the final states with 
T = 5 consist of one (for each type) huge cluster together with a small number of remaining 
agents scattered around. Varying the density v does not radically alter the qualitative structure 
of the final states. We now quantify the aggregation for each value of T, as v varies between 
2% and 33%. 

4.2.1 T = 3: sparse clusters 

The sparsity of the final states with T = 3 is due, in part, to large blocks of the initial checker- 
board configuration that remain unchanged during the evolution. We call this phenomenon the 
super-stability of the checkerboard. Every agent is not just happy, but has one like neighbor 
to "spare". Thus, it takes a large perturbation to make a given agent move and, therefore, only 
agents close to the initially perturbed sites move. Consequently, Schelling required a large den- 
sity of vacant spaces v (33%) to overcome checkerboard super- stability. One can see from the 
left column in Fig.|2]and from Fig.|5]in the supporting materials, that as v decreases, larger and 
larger parts of the initial configuration remain unchanged during the evolution. 
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Figure 3: Statistics of four key measures of aggregation of final states for T = 3 (red triangles), 
T = 4 (green squares), and T = 5 (blue circles) for different v: A The scale of aggregation L; 
B The number of clusters Nc', C The number of agents with eight like nearest neighbors Nq; D 
Normalized perimeter p. 
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We observe that for T = 3, larger values of v result in larger clusters, and thus lead to greater 
global aggregation. The number of clusters in the final states, Nc, decreases as v decreases 
(Fig. |3p) and the dependence is almost cubic (see Fig. |9]in the supporting materials). The 
seclusiveness measure, Nq, is a monotonically decreasing function of v: as v decreases, the 
final state approaches the checkerboard and, naturally, almost all the agents have some contacts 
with other agents. Similarly, the smaller the value of v, the larger the normalized perimeter, p. 

4.2.2 T = 4: compact clusters and mesoscale aggregation 

Increasing T from 3 to 4 eliminates the checkerboard super- stability phenomenon and results 
in strikingly different structures of aggregation in final states. Namely, the final states consist 
of relatively small numbers, that clearly depend on v, of compact clusters (see Fig. [2p -F and 
Fig.|6]in the supporting materials). 

For relatively large v, such final states exhibit mesoscale aggregation and, for small values 
of V, macroscale aggregation. There seems to be no canonical way to separate the two types 
of aggregation. Our criterion to define the transition when the size of the largest cluster, L, 
becomes equal to A^. 

We find two measures that clearly differentiate the global clustering of the final states for 
T = 3 and T = 4. First, the final states have very different perimeters (see Fig.|3p ). Second, 
for T = 3, the clusters are very sparse, while for T = 4, the clusters are compact. A natural 
way to quantify this is to use Nq, whose statistics we present in Fig. [3]C . The measure Nq is a 
monotonically increasing function of v. 

We quantify the differences in the final states for T = 4 with v ranging from = 33% down 
io V = 2%, with three measures: the number of clusters, Nc, the scale of aggregation, L, and 
the seclusiveness, Nq (Fig.|3]). By providing opportunities for increasingly "easier satisfaction," 
one might believe that decreasing v increases the values of Nq. Our study confirms this, and 
the dependence is remarkably linear. Namely, for typical final states with T = 4, t> = 33% 
(Fig. [2^) Nc is relatively high; for T = 4, = 15%, Nc is smaller (Fig. |2^) and the clusters 
on average are bigger; finally, states with T = 4, i; = 2% contain only a few compact clusters 
of either type that stretch across the whole lattice (L = 100). In general, as v decreases, L 
increases almost linearly (see Fig.|3]\). 

Thus for T = 3 and T = 4, the increase in v leads to the opposite effects: they increase and 
decrease the level of global aggregation, respectively. 

4.2.3 T = 5: final states with many unhappy agents 

For small v, the dynamics with T = 5 results in a final state achieved after just a few switches, 
and consists of mostly unhappy agents with no vacant space to where they could move. How- 
ever, a slight modification of the selection algorithm, to allow direct R-B switches when it is not 
possible to switch with a vacant space, results in significant global aggregation and drastically 
reduces the number of unhappy agents, although not eliminating them entirely. The presence 
of unhappy agents in the final states is a new phenomenon, which we do not observe in simula- 
tions for T = 3 (while such configurations theoretically exist, they are extremely unlikely) and 
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is much less pronounced for T = 4 (see Fig.[8]\ in the supporting materials). 

A typical T = 5 final state with modified selection algorithm consists of one big cluster for 
each kind of agent and the rest of the agents are unhappy and scattered around. The smaller that 
the value of v is, the larger the two main clusters are and the fewer unhappy agents there are. 

The globally aggregated final states (small values of v) with T = 5 (with modified selection) 
and T = 4 (with Schelling selection) appear similar in terms of the number of large clusters 
and the scale of aggregation, L (Figs. |2p and G). However, there is a large difference in their 
perimeter p: it is much smaller for T = 5 (see Fig. [3p). 

5 Concluding Remarks 

In this paper we devised new measures to quantify the aggregation in the Schelling segregation 
model and studied their dependence on city size, disparate neighbor comfortability threshold, 
and population density. We identified distinct scales of global aggregation, and showed that 
the striking global aggregation Schelling observed for the 8x8 lattice is strictly a small lattice 
phenomenon. We also discovered several remarkable scaling laws for the aggregation measures 
as functions of population density. 

6 Supporting Material 
6.1 Schelling's Simulations 

Schelling considered the cases = 8, T = 3, and v = 33%. For T = 3 or 4, and V = 0, a 
"checkerboard" configuration of B's and R's (imagine placing B's on the red squares and R's 
on the black squares of an actual checkerboard) is a final state, since all agents have four like 
nearest neighbors. 

To generate his initial configurations, Schelling begins with a checkerboard configuration 
without periodic boundary conditions and randomly removes approximately one third of the B's 
and R's, keeping equal numbers of both agents (9). We refer to the result as a deleted checker- 
board configuration. Removing these agents makes some of the remaining agents unhappy and 
drives the evolution. Several authors have observed that removing such a large percentage of 
agents is unnatural, but it is crucial to attain aggregation in Schelling's model. Removing fewer 
agents results is a final configuration close to the initial configuration. 

Finally, to construct his initial configurations, Schelling modifies the deleted checkerboard 
by randomly adding a total of 5 B's and 5 R's in vacant spaces. Such an initial configuration is 
assumed to be a proxy for a nearly integrated city. 

A typical run of Schelling's original model system (with periodic boundary conditions) is 
presented in Fig. |4j initial state (Fig. |4]\) and final state (Fig. |4p). Schelling performed many 
simulations by hand using an actual checkerboard, and observed that the final states presented 
a significant degree of global aggregation. He equated the global aggregation with segregation 
of a city. 
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Figure 4: Schelling's original model with N = 8: A: initial state and B: final state. 

6.2 Structure of final states for T = 3 

In Fig. [5] we present typical final states for T = 3 and different values of v. The pictures and our 
statistics show that as the value of v decreases, the number of clusters decreases (Fig.|3p), the 
final states retain larger and larger parts of the original checkerboard configuration. One can see 
from Fig.[5]that, unlike the final states for Schelling's original model with = 8 (see Fig.|4p) 
that consists of just one or two separate domains of R and B agents, the final states for N = 100 
and V = 33% contain many clusters. The striking qualitative difference is also quantified by the 
relatively large values of the normalized perimeter, p (see Fig.[3p) and by the large differences 
in the values of Nq/N"^ between the states presented in Fig. [3^ and |4^. 



6.3 Structure of final states for T = 4 

In Fig. 6 we present typical final states for T = 4 and different values of v. The pictures and 
our statistics show that as the value of v decreases, the number of clusters decreases (Fig.|3p) 
and the size of the compact clusters (scale of aggregation) increases (Fig.[3]\). 



6.4 Structure of final states for T = 5 

In Fig. [7] we present typical final states for T = 5 and different values of v. Observe that most 
of the agents away from the big clusters are unhappy (see Fig.[8]for the statistics of the unhappy 
agents). The pictures and our statistics show that as the value of v decreases, the number of 
unhappy agents in final states v decreases (Fig. [8]\), and the size of a single (for each type) 
major cluster increases (see Fig. [8p). Another clear indication of the growth of the main cluster 
is the increase of N^, presented in Fig. |3p. 

The average number of unhappy agents in final states for different values of T and v is 
presented in Fig. [8]\ . It is remarkable that the average number of unhappy agents is almost a 
linear function of v, between v = 10% (where they constitute approximately 10%) and v = 30% 
(where they constitute approximately 33%, in other words, almost every agent). While the 
existence of unhappy agents in the final state does not significantly increase the perimeter p of 
the final states, it greatly inflates the total number of clusters Nc- 
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Figure 6: Characteristic final states for T = 4 for different v: A: v = 2%, B: v = 6%, 
V = 10%, D: ^; = 15%, E: v = 20%, F: v = 24%, G: v = 28%, U v = 33%. 
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Figure 8: Statistics of the final states with T 

in final states for T = 3 (red triangles), T = 4 (green squares), and T 
average number of the agents in the two big clusters. 



5: A: The average number of unhappy agents 

5 (blue circles); B: The 



10= 




JQ-i JQ-I ](,0 

log(v) 

Figure 9: Statistics of the number of clusters in the final states for T = 3 (red triangles), T = 4 
(green squares) for different f on a log — log scale. 

6.5 Statistics of the number of clusters 

In Fig. [9] we replot the data shown in Fig.[3p on a log scale (to illustrate power-law statistics), 
showing the statistics of the number of clusters in the final states for T = 3 and T = 4, for 
different v. The lines indicate power laws. The values for the slopes are 0.89 for T = 4 (the 
value 1 is well within the error bars) and 2.76 for T = 3 (the value 3 is well within the error 
bars). 



6.6 Number of steps in the evolution 

The average number of switches required for the final state to be achieved are (3596, 5192, 
5573, 4422) respectively, for T = 3, i; = 33%; T = A,v = 33%; T = A,v = 2%; and 
T = 5,f = 2%. The switches for T = 5,v = 0.02 that included both R and B agents are 
counted twice. The most striking feature is that it takes significantly fewer switches to achieve 
the final state for T = 5 than for T = 4. The distribution of the number of jumps for different 



agents is presented in Fig. 10 
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Figure 10: Statistics of the number of agent jumps. 

6.7 Final states with N = 50 and N = 200 

To illustrate the dependence of the final states of A^, we performed the simulations for = 50 
and for = 200. Characteristic final states are presented in Fig.fTTland Fig.[T2]for = 50 and 
A^ = 200, respectively. One can see that both figures are qualitatively very similar to Fig. |2] 



6.8 Measures of Aggregation 

Schelling considered two quantities to measure the aggregation of a state: 

1. The ratio of unlike to like neighbors is called the [u/Z] -measure. For a site on the final lattice 

with coordinates we define: 

[u/^kj = — — 



where Pij, qij, and Wij are the number of like, unlike, and vacant neighbors of the agent located 
at (i, j), respectively. We define the sparsity ([«//]) of a cluster by averaging the [m//] -measure 
over the given cluster. 

2. The number of agents that have neighbors only of the same kind, A'o (note, that this definition 
excludes the vacant spaces as well). The abundance of such agents indicate the presence of 
large, "solid" clusters. This quantity is the most useful in quantifying between the states with 
T = 3 and T = 4. 

Besides those two characteristics, we introduce several more. 

3. The total number of clusters, A',^. With 8-point neighborhoods, clusters may be intermingled. 
For example, a checkerboard configuration has just 1 + 1 = 2 clusters. The quantity Nc is 
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Figure 12: Characteristic final states for = 200 and different values of T and v: A: T = 3, 

V = 2%, B:T = 3,v = 15%, C:T = 3,v = 33%, D:T = A,v = 2%, E:T = A,v = 15%, F: 
T = 4:,v = 33%, G:T = 5,v = 2%, H: T = 5, i; = 15%, l:T = 5,v = 33%. 
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the most useful for configurations consisting of compact clusters of a similar size. To study 
configurations such as the final states for T — 5, a more useful quantity is the number of 
clusters that have greater then, say, M^aj./10 agents, where M^^ax is the number of agents in 
the largest cluster of a given kind. 

4. Total perimeter (normalized by the number of agents). For 7^ 0, we follow the calculation 
of surface tension in the physics of foams, and write 

P 1 ^ 

Thus defined, P plays the role of a Lyapunov function: it can be shown that in the process of 
evolution every legal switch makes P smaller. 

Suppose we switch an R and a V. Before the switch, suppose R had Bi, Ri, and Vi, of B, R, 
and V neighbors, respectively. Similarly, the numbers for the V agent are D2, R2, and V2. Then 
the value of P before and after the switch are: 

Pinitial — 2-Bl + Vi + B2 + R2; P final — 2i?2 + V2 + -Bi + -Ri- 

Thus, 

P final — Pinitial — B2 + V2 + R\ — {Bi -\- Vi -\- R2) ■ 

Taking into account that Bi -\-Vi -\- Ri — B2 -\-V2 -\- R2 — 8, we arrive at 

P final — Pinitial = 2 {Ri — R2) < 0. 

Similarly, if the switch between R and B, we have 

Pfinal — Pinitial = 2 (-R2 — Rl) + 2 — B2) < 0. 

5. We define the diameter of a cluster, L, to be the side length of the smallest square that covers 
the cluster. The diameter of a cluster can be easily computed as the larger between the number 
of rows that contain an agent belonging to the cluster and the number of columns that contain 
an agent belonging to the cluster. For configurations consisting of mostly compact clusters, the 
maximum diameter, L = max(Lj), defines the scale of aggregation. 

References and Notes 

[1] 1. Benenson, E. Or, E. Hatna, I. Omer, Residential Distribution in the City Reexamined, 
in press. 

[2] E. Bruch, R. Mare, Neighborhood Choice and Neighborhood Change, AIS 112, 667-706, 
2006. 



18 



[3] M. PoUicott and H. Weiss, The Dynamics of Schelling-Type Segregation Models and a 
Nonlinear Graph Laplacian Variational Problem, Advances in Applied Mathematics 27, 
17-40 (2001). 

[4] lhttp://www.realestatea gent.com/glossary/real-estate-glossary-show-term-1699- 
occupancy-ratio .html , 

[5] B. Simon, The Statistical Mechanics of Lattice Gases, Vol 1, PUP, 1993. 

[6] T. Schelling, Models of Segregation, American Economic Review 59, 488-493, 1969. 

[7] T. Schelling, On the Ecology of Micromotives, The Public Interest 25, 61-98, 1971. 

[8] T. Schelling, Dynamic Models of Segregation, Journal of Mathematical Sociology 1, 143- 
186, 1971. 

[9] T. Schelling, Micromotives and Macrobehavior (W.W. Norton, New York, ed.l, 1978), 
pp. 137- 166 

[10] D. Vinkovic and A. Kirman, A Physical Analogue of the Schelling Model, PNAS 103, 
19261-19265, 2007 

[11] S. Gerhold, L. Glebsky, C. Schneider, H. Weiss, and B. Zimmermann, Limit States for 
One-dimensional Schelling Segregation Models , Commin Nonlinear Science and Numer- 
ical Simulations, to appear. 

[12] A. Singh, D. Vainchtein, and H. Weiss, Schelling' s Segregation Model with Minorities, In 
preparation. 

[13] W. A. V. Clark, Residential Preferences and Neighborhood Racial Segregation: A Test of 
the Schelling Segregation Model, Demography, 28, 1-19, 1991. doi: 10.2307/2061333. 

[14] R. Sander, D. Schreiber, J. Doherty, Empirically Testing a Computational Model: Archi- 
tectures and Institutions, edited by D. Sallach and T. Wolsko, Chicago, IL: University of 
Chicago; ANL/DIS/TM-60, Argonne National Laboratory, Argonne, IL, 2000. 



19 



